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Abstract. We introduce a model for active transport on inhomogeneous networks em- 
bedded in a diffusive environment and investigate the formation of particle clusters. 
In the presence of a hard-core interaction, cluster sizes exhibit an algebraically decay- 
ing distribution in a large parameter regime, indicating the existence of clusters on all 
scales. The results are compared with a diffusion limited aggregation model and active 
transport on a regular network. For both models we observe aggregation of particles 
to clusters which are characterized by a finite size-scale if the relevant time-scales and 
particle densities are considered. 

PACS. 87.16.Uv Active transport processes - 45.70.Vn Granular models of complex 

systems 



1 Introduction 

Active transport processes play an important role 
for the functionality of social and biological sys- 
tems. Examples are road traffic or active intra- 
cellular transport of vesicles and organelles by 
motor proteins that perform directed movement 
along microtubules or actin filaments. In recent 
years stochastic systems of self-driven particles 
have been applied to model intracellular trans- 
port of motor proteins in a number of works 
|l|2|3|4|5|6)7j . While the microtubules usually ar- 
range in an ordered pattern (e.g. a radial struc- 
ture in most mammalian cells, though longitu- 
dinal in neuronal axons), actin filaments often 
form randomly structured undirected networks. 
Therefore the investigation of transport on net- 
works arises to be an interesting object of re- 
search. 

Active transport on an undirected but reg- 
ular network has been investigated by Klumpp 
et al. [8]. They studied the dynamic properties 
of diffusive non-interacting particles on a lat- 
tice with an embedded regular square network 



consisting of active stripes where particles per- 
form biased motion. They showed that though 
movement of particles remains globally diffusive 
on long time-scales, the diffusion constant is en- 
hanced by the presence of the network. 

In transport systems considering steric exclu- 
sion interaction between particles, aggregation is 
a common phenomenon, manifesting in the for- 
mation of jams. Jams can form in one dimen- 
sional systems with single tracks due to bound- 
ary conditions {boundary induced phase transi- 
tions 0), induced by defects [10111112113114] or 
they emerge spontaneously due to stochastic slow 
down of vehicles in highway traffic |15|16j . In two 
dimensional regular street networks, mutual in- 
terference of vehicles at intersections lead to jam- 
ming [IT]- 

Studies of transport on inhomogeneous topo- 
logical networks (graphs with nodes and edges, 
no distances) revealed an interesting phenomenol- 
ogy. E.g. non-interacting particles performing a 
random walk exhibit an inhomogeneous density 
distribution on the nodes [18] , while inclusion of 
an attractive zero range interaction even allows 
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the particles to aggregate and form a conden- 
sate, corresponding to nodes containing a finite 
fraction of particles [19 . These results show that 
the structure of a transport network strongly in- 
fluences transport properties. In order to model 
active transport on act in filament networks, it is 
therefore necessary to consider realistic network 
structures. 

For many biological processes, concentration 
gradients are crucial. One example is the ag- 
gregation of proteins inside the cell or in the 
cell membrane. Clusters of aggregated proteins 
can be observed and characterized experimen- 
tally for example by fluorescence microscopy. In 
some cases these clusters are essential for cell 
functionality but they can also lead to dysfunc- 
tions or even apoptosis. In yeast cell membranes 
for example one observes the aggregation of Erd2p- 
receptors which can promote the internalisation 
of toxins [20^. The existence of a quasi two- 
dimensional irregular actin filament network be- 
neath the membrane |21|22j suggests a jamming 
mechanism of vesicles prior exocytosis, resulting 
in clusters of receptors on the membrane surface. 
The limits of resolution in optical microscopy 
[23] do not allow the localization of single recep- 
tors. By contrast the size of larger particle ag- 
gregates can in principle be given with relatively 
high precision. Therefore it is useful to relate the 
cluster size distribution with microscopic trans- 
port mechanisms by means of theoretical mod- 
elling. 

In this work we propose a model for active 
transport of extended hard-core particles (corre- 
sponding to vesicles) on a two-dimensional ran- 
domly disordered network embedded in a diffu- 
sive environment. The model is motivated by in- 
tracellular transport on submembranal networks, 
we therefore adapt the model parameters to this 
reference system. We check particle configura- 
tions in order to identify the formation of clus- 
ters and investigate cluster size distributions. The 
results are compared with a regular network in 
diffusive environment and a diffusive system with- 
out network where attractive particle-particle in- 
teractions promote cluster formation similiar to 
a van der Waals gas. The main focus will be on 
robust properties of clusters that serve as crite- 



ria to discriminate between different microscopic 
dynamics. 



2 Network models 

In the following, we introduce stochastic mod- 
els in order to study the influence of the net- 
work structure on dynamical properties. While 
the presence of a network in a diffusive envi- 
ronment without interaction between particles 
merely leads to the enhancement of the diffusion 
process [8 , we are rather interested in collec- 
tive patterns emerging due to interactions. Our 
simulations use stochastic dynamics in order to 
integrate the many-particle Master equation. At 
each time step particles are randomly chosen and 
updated (random sequential update) according 
to the rules given in table [l] applying periodic 
boundary conditions. Time steps are normalized 
so that on average each free particle performs 
one diffusive step per time step At. 



2.1 Regular networks 

As a the first example for active transport on 
networks a discrete lattice gas model similiar 
to the model investigated in [8 is considered. 
N X N sites are arranged in a square lattice of 
edge length Na where a is the lattice spacing. 
Each site can either be empty or occupied by 
at most one particle. We distinguish the particle 
states attached {A) and detached (D). The system 
contains stripes of active sites that constitute a 
regular square transport network where attached 
particles perform a directed motion. Detached 
particles always move diffusively. The orienta- 
tion of stripes was chosen randomly with equal 
probability. Attached Particles on active stripes 
can unbind becoming detached, while detached 
particles can bind to stripes becoming attached. 
Steps that would result in double occupation of 
a site are prohibited if not state differently. 

Compared to the dynamics of non-interacting 
self-driven particles qualitatively new features 
arise due to the steric particle-particle interac- 
tions at intersections of the network. Here we 
introduce an additional parameter, the blocking 
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Fig. 1. Illustration of the dynamics in a regular net- 
work. Dark grey discs are particles diffusing freely 
with rate cod- Black discs represent attached par- 
ticles that can only step in the preferred direction 
of the active stripe they occupy (bold arrows) with 
rate p. On filament sites (light grey), particles can in- 
terchange between attached and detached state with 
rates coa and cOd, respectively. Crossed arrows denote 
steps that are inhibited due to the exclusion princi- 
ple. 



tides can perform directed paths along these fil- 
aments. 



2.2.1 General properties of the model 

The main components of our model are filaments 
and particles interacting via a spherical hard- 
core potential represented by a disc of radius r^. 
This hard core potential is implemented by can- 
celling any steps that would result in an overlap 
of discs. Filaments are represented by linear se- 
quences of nodes with a distance of dn between 
nodes. They are directed with a (-)-end at which 
they can shrink and a (+)-end at which new 
nodes can be generated to elongate the filament. 
Particles can attach to nodes that are within a 
distance less than d^ and perform steps to adja- 
cent nodes in the (+)-direction of the filament. 

2.2.2 Dynamics of filaments and particles 



probability: If at least two particles are at sites 
adjacent to an intersection site, each particle may 
only access the intersection site with the proba- 
bility 1 — b (cf. figure [T]). Particles on intersection 
sites retain their moving direction. 

The explicit rules for the particle dynamics 
are displayed in table [T] and illustrated in figure 
[l] We have chosen the default parameter values 
analogous to ^ lod = P = 0.5, cOd = 0.02 and 
system size N = 200. In [8 the attachment rate 
is equal to one, which corresponds to an effective 
attachment rate cOa = 0.25 if a particle is on 
an adjacent non- active sit^ To be consistent 
with the subsequent continuous space model, we 
choose 6 = 1. 



2.2 Inhomogeneous networks 

Generalising to continuous space and allowing 
for arbitrary directions and lengths of active stripes 
we present a continuous model with randomly 
generated linear filaments where hard-core par- 



^ In contrast to 8 , we allow for crossing of active 
stripes by diffusion. 



The filament network is generated by stochas- 
tic dynamics, updating network configurations 
by the processes commented in table [2] and illus- 
trated in fig. |3] Since the model is motivated by 
the dynamics of actin filaments some additional 
parameters originating from the biological pro- 
cesses in cells are included. 

The quantity Parp introduced in table [2]rep- 
resents the density of free ARP2/ 3- complexes 
that serve as nucleation and branching seeds for 
filaments, while the actin density pact corresponds 
to the density of free nodes (actin monomers) 
constituting the filaments. Their initial values 
are p^Rp Pact which corresponds to the case 
if all monomers are dissociated. The densities 
decrease with the growing filament network as 
shown in fig. |4] After 5000 steps a stationary 
actin density is reached. We assume that also 
the structure of the network is stationary and 
in qualitative agreement with a real actin net- 
work. We therefore stop network dynamics at 
this point. 

After construction of the network particles 
obeying the exclusion principle are fed into the 
system at random positions. As mentioned above, 
the particle positions are updated following a 
random sequential update scheme, whereby the 
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Process 



Particle state (s) Description 



Probability 



Diffusion 
Forward Step 

Attachment 
Detachment 
Blocking 



D 

A 

D^A 
A^D 
D 



Detached particles move to sites randomly chosen from cji) = 1 
the four neighbors 

Attached particles move to the next site in forward di- p = 0.5 
rection of filament 

Detached particles on filament sites becomes bound ooa = 0.25 

Attached particles become detached ood = 0.02 

Forward movement of particles adjacent to intersection 6=1 
sites is inhibited if other particles occupy sites adjacent 
to intersection 



Table 1. Brief prescription of the dynamic processes in the square lattice model. Column 2 displays 
the particle states "attached" (A) and "detached" (D). The right column displays the probability that the 
respective process occurs within one time step. Numerical values given in the right column are default values 
which are used if not stated else and are chosen to fit the ones in [8]. 



Process 



Description 



Probability 



Nucleation 



Branching 



Growth 
Shrinking 

Uncapping 



Initialization of filaments with arbitrary direction at an 
arbitrary point in the system. The (-)-end receives a cap 
inhibiting shrinking. 

New filaments are initialized at an existing one (not nec- 
essarily the (+)-end; angle between parent filament and 
branch=70^24 . 

New nodes are generated at the (+)-ends of filaments 
Nodes are removed at the (-)-end of filaments if the end 
is not capped. 
Caps are removed. 

Table 2. Dynamics of the filament network. 



OJn pact pARP 



Pact PARP 



(^gpact 
00s 



particle-particle as well as the interactions be- 
tween particles and the generated static network 
are considered . Like in the regular network model 
particles can freely diffuse in the 'detached' state 
and perform directed movement in the 'attached' 
state. The rules of the particle dynamics are pre- 
scribed in table |3] and illustrated in fig. [2] 

Although we do not consider a particular bi- 
ological system, we choose parameters to fit the 
typical order of magnitude in real vesicular trans- 
port. If not stated differently, we will use default 
parameters displayed in table [4] for our simu- 
lations. The referenced works used experimen- 
tal and modeling techniques to obtain the data 
given in the third column. For particle dynamics, 
we choose the parameters to be consistent with 



the discrete model introduced in the last section 
relying on the model in 

3 Results 

3.1 Characterisation of Clusters 

As already mentioned in the introduction we aim 
to relate the microscopic particle dynamics to 
the size distribution of their aggregates. In this 
section we discuss the definition of clusters for 
the different model systems. 

Clusters are groups of particles that are con- 
nected by overlapping neighborhoods. We there- 
fore introduce the X— neighborhood of a parti- 
cle representing a disc of radius A around the 
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Process Particle state (s) Description Probability 



Diffusion 


D 


Detached particles move in a random direction. Step 
widths are uniformly distributed between and 21 d 




Step 


A 


Attached particles move to adjacent node in (+)- 
direction. 


P 


Attachment 


D^A 


Particles bind to nodes if their distance is less than db, 
becoming 'attached' 




Detachment 


A^D 


Particles detach 





Table 3. Particle dynamics. A= 'attached'; D= 'detached' 



Parameter name 


Reference 


Reference Value 


model parameters 


Filament dynamics: 








nucleation rate Un 


M 


8.7 X IQ-^pM-^s-^ 


1.0 X 10-^lu-^tu- 


growth rate oOg 


M 




0.25 lu^ tu-^ 


shrink rate ujg 


[26] 


4.2s"^ 


0.075 tu-^ 


branch rate Ub 


[25] 


5.4 X 10"^ /xM-^s-^ 


0.0001 lu^ tu-^ 


uncap rate Uu 


[25] 


0.0018 s-^ 


O.OOOltu-^ 


actin density pact 




free actin: 0.1 — IpM 


2lu-^ ^ 


ARP2/3 density p%p 


^ 


O.lpM 


0.1 lu-^ 


Particle dynamics: 








particle radius rp 


m 


42.5nm (average) 


0.5 lu 


binding distance db 


m 


1 site (50nm) 


0.5 lu 


node distance dn 


m 


36nm 


0.36 lu 


attachment ooa 


m 


1/4 of diffusive steps 


0.25 tu-^ 


detachment Ud 


m 


0.8s-^ 


0.02 tu~^ 


diffusive step Id 


m 


1 per time step 


0.5 lu 


step rate p 


El 


20s~"^ ^ V = 1pm/ s 


0.75 tu-^ 


particle density p 




10-60 vesicles in bud (^0.75/im radius) 


0.04lu-^ 



Table 4. Default parameters of the model which are biologically motivated by transport of vesicles 
by myosin on actin filaments. The referenced values are either based on experimental data or existing 
models for intracellular transport and filament dynamics 28 . Model parameters are chosen to be in 
the order of magnitude of referenced values, fitted to time and space scale of the simulations. Length 
scale: llu = lOOnm = 2rp IpM = 0.6 Time scale: Itu — At — 0.025s. By default we consider 
square systems of system length L = 200 lu. 

^ As shown in figure ^ the amount of overall actin is about ten- fold larger than free actin. Thus we 
chose this factor to extrapolate the experimantal date to p^act 



center of the particle. A cluster is defined as a 
set of particles included in a connected area of 
A— neighborhoods (cf. fig.|5|. If continuous space 
variables are used, there exists no natural scale 
which identifies two particles as neighbors. We 
therefore have to specify the value of A. In or- 
der to extract relevant results, we fix A so that 
qualitative results are robust on deviations. If 



not stated differently we choose A = 2r^, which 



turns out to meet this condition (cf. fig 19). 

In lattice models, static particle clusters are 
usually considered as as connected sets of ad- 
jacent particles. However this definition is not 
appropriate in this context since clusters move 
by propagation of vacancies. Therefore we con- 
sider particles separated by a single vacancy as 
belonging to the same cluster. 
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Fig. 2. Illustration of the particle dynamics in an in- 
homogeneous network. Dark grey discs are free par- 
ticles, black discs represent particles attached to fila- 
ments stepping to adjacent nodes (distance dn) with 
rate p. Particles can attach to filaments with rate Ua 
if they are within the binding distance db and de- 
tach with rate ood- Overlapping is inhibited due to 
exclusion. 




Fig. 5. Illustration of particle clusters. Black 
discs represent particles, while grey discs are the A- 
neighborhoods of each particle. Connected grey ar- 
eas are clusters; the size of a given cluster is the 
number of particles on it. 



(-)-End (+)-End 



Branching 



Fig. 3. Illustration of the filament dynamics. Fila- 
ments are implemented as sequences of nodes (small 
dots, corresponding to actin monomers) separated 
by a distance dn (short bars) . Filament segments are 
polarized, with a (+)-end where nodes are created to 
elongate, and a (-)-end where nodes dissociate caus- 
ing shrinking. 
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time steps 

Fig. 4. Density of free actin 
on time. After 5000 time steps : 
reached. 



10000 

pact in dependence 
I stationary state is 



Our main interest is in ensemble and time 
averages of cluster size distributions (CD) and 
their asymptotic behaviour. CDs display the rel- 
ative frequency of cluster sizes emerging in the 
system. If not stated differently we averaged over 
50000 time steps within individual runs, evaluat- 
ing cluster distributions in distances of 500 time 
steps, taking an ensemble of 100 samples. 

Clustering also occurs for random particle 
configurations. In fig. |6] and [7] we displayed clus- 
ter size distributions of random configurations in 
discrete and continuous space for different par- 
ticle densities. Here the density p is the parti- 
cle number per area unit which corresponds to 
(2rp)^ in continuous space and one site in the 
lattice model. If densities are not too large, the 
formation of large clusters is impeded resulting 
in an exponentially decaying cluster size distri- 
bution. For high densities one observes a small 
peak at the right end. At these densities clusters 
spanning the whole system emerge. In order to 
rule out these kinds of random clustering we will 
only consider densities below the regime of span- 
ning clusters at relevant scales A. In this work we 
are interested in cluster formation mechanisms 
beyond random clustering. 

In the following, we will focus on particle con- 
figurations and cluster-size distributions in sev- 
eral transport models. 
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Fig. 6. Cluster size distributions of random particle 
configurations in dependence on the particle density. 
For low densities the CD decays fast with a short 
size-scale. For large densities clusters on large size- 
scales and even such that span the whole system 
emerge (not visible in figure since on too large size- 
scale) . 



(a) continuous space 

ir 
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(b) discrete space 
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Fig. 7. Cluster size distributions of random particle 
configurations for different coarse graining scales. A 
is the radius of the environment as defined in section 



3.1 while d represents the maximum distance (num- 
ber of vacancies) allowed between two particles con- 
necting a cluster. One observes that for large coarse 
graining scales clusters spanning the whole system 
emerge. 



3.2 Aggregation without network 

A possible scenario for clustering of receptors on 
the cell surface is aggregation due to an attrac- 
tive interaction between proteins diffusing in the 
cell membrane. This process can be formulated 
as an equilibrium model consisting of diffusing 
hard-core particles (radius Vp ~ lOnm; within 
the size-scale of membrane proteins) interacting 
via an attractive potential. We app ly the par- 
ticle dynamics discussed in sec. |2.2| but do not 
consider filaments. In addition we introduce a 
particle-particle interaction between the parti- 
cles realized by a square well potential of the 
form 



l/(x - xO 



-Vo for |x-x^| < 







for 



> dv 



(1) 



where x, x' are particle positions. This potential 
can be implemented using a Metropolis accep- 
tance probability p = min(e~^'^^'^^'^+i^~^'^^'^^\ 1) 
for a step from Xn to x^+i (n denotes the time in- 
dex) . In the following we use dimensionless quan- 
tities and put P = 1. The default parameters are 
Vo = 2,dv = 3.5rp and particle density p = 
0.04. Assuming a diffusion constant for mem- 
brane proteins D ?^0.0025/im^/s [31 we choose 
a time step corresponding to At = 0.02 seconds 
so that one diffusive step of length Vp = lOnm is 
performed per time step At. 

In fig. [8] typical particle configurations at sev- 
eral runtimes are displayed, while in fig. |9] en- 
semble averages of cluster size distributions are 
shown. Initial clustering already occurs on a rather 
small time-scale. Regarding the particle config- 
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urations we see that the number of clusters de- 
creases with increasing runtime while the aver- 
age size of remaining clusters increases. This is 
due to diffusion and merging of existing clus- 
ters after long times. Movement of large clus- 
ters is strongly supressed, so that merging oc- 
curs quite slowly. The coarsening process can 
also be observed in the cluster size distribution. 
We observe a characteristic scale for larger clus- 
ters, manifesting in the emergence of a maxi- 
mum, indicating a characteristic scale for cluster 
sizes. The dominant clusters always are within 
the same size-scale which increases with time. 

Since the cell membrane changes its struc- 
ture continuously, patterns arising at time-scales 
corresponding to a finite fraction of a cell cycle 
cannot be assumed to be in a stationary state. 
Computing time averages we therefore focus on 
intermediate times and fix the averaging interval 
starting at 20000 time steps (corresponding to 
r^l minutes in real time) after random initialisa- 
tion of particles and ending at 30000 time steps. 
The time interval lies in the transient regime for 
default parameters. Within this interval we com- 
puted cluster size distributions (time and ensem- 
ble averages, 200 samples) for different param- 



eter regimes and displayed them in figure [10 
One observes that for weak interaction no sig- 
nificant clustering occurs manifesting in an ex- 
ponentially decaying CD, while for strong inter- 
action Vq^ including the default parameters, a 
maximum emerges hallmarking the formation of 
clusters. 



3.3 Directed transport on regular networks 

In this section we examine features of particle 
configurations and cluster distributions in the 
model introduced in section |2.1[ i.e. a regular 
network of active stripes. As in the last section 
we start time averaging after ts = 20000 time 
steps. We carefully checked that a stationary state 
has been reached at this point, (cf. fig. 12). As 
time averaging interval we choose 50000 time 
steps. In fig. [TT] particle configurations for mod- 
erate and high densities are displayed. For parti- 
cle density p = 0.04 one observes small L-shaped 
clusters centering at intersections. For higher den- 
sities it appears that clusters that are becoming 
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Fig. 8. Configurations of particles (black discs 
= particle neighborhoods with radius A = 2rp), 
exhibiting a mutual attractive interaction. Snap- 
shots at different times for Vb = 2,dv = 3.5, p = 
0.04, L = 200, 1 timestep=0.025sec. One observes 
that already at small times clusters form and for 
long run times they the number of clusters decreases, 
while the average size of cluster increases. 
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Fig. 9. Cluster distributions in dependence on the 
runtime in a system without network but attrac- 
tive square well interaction potential. Parameters are 
Vb = 3, = 3.5, p = 0.04, L = 200 r^, average over 
200 runs. A maximum establishes, that moves slowly 
towards larger scales. 



larger and merge with each other to form large 
mesh-shaped clusters (cf. fig[TTJb)). However, in 
this case clusters are hardly distinguishable and 
not well separated which results in sensitive de- 
pendence on the coarse graining scale (cf. fig. 
14). In fig. 13 we plotted the cluster size distri- 
butions averaged over time and 100 individual 
runs. Examining the cluster size distributions in 
fig. [Tsj one observes an exponential decay for 
densities which are biologically relevant (cf. see 
also the configuration in figure 11 ^a)). Therefore 
only small clusters emerge while large clusters 
are exponentially supressed. 

For large densities (> 0.1) the decay of the 
cluster size distribution becomes algebraic, indi- 
cating that clusters on all size-scales exist. These 
large clusters correspond to the ones generated 
by merged small clusters as displayed in fig.jlTJb). 



3.4 Inhomogeneous networks 

The filament growth dynamics described in sec- 
tion |2.2| generate a network where single fila- 
ments have random length and direction. In or- 
der to keep dynamics simple but retaining the 
crucial features of disordered networks, we here 
neglect branching and the dynamics of ARP2/3 
to obtain a network of uncorrelated filament ori- 
entations. The other processes are required to 
obtain a disordered stationary network configu- 
ration. 




20 40 60 

cluster size 




Fig. 10. Plots of cluster distributions in the at- 
tractive particles model at intermediate times (t — 
20000 — 30000 time steps) in dependence on the po- 
tential depth Vb (a), the potential width dy (b) and 
particle density p (c). One observes the transition 
from an exponential decay (non clustering phase) to 
the formation of a maximum, corresponding to clus- 
ters at this size-scale (condensation) . Default param- 
eters are given in the text. 
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(b) y9 = 0.15 




Fig. 11. Particle configurations in a regular square 
network with system length N — 100 sites. The black 
discs represent particle neighborhoods with radius 
1.1 sites so that discs of particles with one vacancy or 
less in between overlap. One observes the formation 
of small L-shaped clusters at intersection points for 
moderate densities (a). For large densities clusters 
merge forming cluster meshes on all size-scales (b). 



3.4.1 Dynamics of interacting particles 



We obtained particle configurations and cluster 
size distributions applying steric interactions, 
which are shown in figures [T5][l9l The time evolu- 
tion of the cluster size distribution (fig.[l7|) shows 
that a stationary state is reached after 10000 
time steps. Starting time averaging after 20000 
time steps (averaging interval=50000 time steps) 
therefore captures the steady state dynamics. 




number of time steps=100 
number of time steps=1000 
number of time steps= 10000 



^^'^-^ 10 20 30 

cluster size 

Fig. 12. Cluster size distributions in a regular net- 
work for different runtimes (given in time steps). For 
a given runtime, we chose the last 100 steps to per- 
form the measurement, taking 100 samples. The CD 
does not change after 1000 time steps, indicating 
that the system is in a stationary state. 



The configuration for p = 0.04 (fig.fTsl shows 
that well separated compact clusters exhibiting 
different sizes emerge (see also scaling in fig. 19), 
while at p = 0.008 no clusters are observed. In a 
large parameter regime including the biological 
relevant default parameters (table[4|, the asymp- 
totic decay of the CD is algebraic in contrast to 
the predominant exponential behaviour on a reg- 
ular network. For a given combination of network 
structure and particle dynamics the properties of 
the CD are determined by the particle density, 
while varying other parameters does not lead to 
qualitative changes except in extreme regimes. 
In fig. [20] cluster size distributions of a regular 

One 



20] cluster size distributions of a ree 
and inhomogeneous network are comparecrl 
observes that clustering is significantly enhanced 
in the inhomogeneous network. 

Due to the finite number of particles there is 
a cut-off at the upper end (e.g. in fig. 19). Fig. 



|l8[a) shows that for increasing system size L the 
cut off regime tends to larger values, indicating 
that this indeed is a finite size effect and asymp- 
totically algebraic behaviour prevails in the ther- 
modynamic limit. Though the exponent 7 of the 
algebraic decay varies for different particle densi- 
ties, the algebraic form is a robust feature. This 
indicates that in the thermodynamic limit clus- 
ters on all size-scales exist. In contrast to reg- 



^ Differences in effective rates due to the different 
spatial character of the system (discrete and contin- 
uous) are not significant since dependence on these 
parameters is weak. 
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p=0.08 
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exponential fit ~ exp(-x/Xj, 



20 40 60 

cluster size 




Fig. 13. Cluster distributions in a lattice gas model 
with exclusion interaction and a regular square net- 
work of active stripes in dependence on particle den- 
sity yo, logarithmic plot (a), double logarithmic plot 
(b) blocking rate b (c),and mesh size / (d). One ob- 
serves that cluster size distributions decay exponen- 
tially for moderate densities resulting in a finite size- 
scale of clusters. For very large densities decay is al- 
gebraic indicating the emergence of clusters on all 
size-scales (see also fig. 11). Default parameters: see 
table H 




10 100 
cluster size 

Fig. 14. Cluster size distributions in the regular 
network obtained by using different definitions of 
the coarse graining scales, d is the distance allowed 
between two particles to connect a cluster. One ob- 
serves that while variations are small for moderate 
density, there is a significant influence on the scale 
for larger densities, indicating that clusters are not 
well separated. 



(a) p = 0.04 




Fig. 15. Particle configurations (black discs = 
neighborhoods with A = 2rp) for default parameters 
{p = 0.04) and system size L = 200 r^. One observes 
big and small clusters. This general picture is pre- 
dominant for a large parameter regime and moderate 
densities. 



ular networks, scale free clustering occurs even 
for moderate densities {p ^ 0.012) exhibiting a 
pattern of well separated clusters. 

The long range behaviour of the CD vanishes 
if crucial properties for cluster formation are ab- 
sent, which can be seen in the particle configu- 
rations in fig. [16] E.g. for small densities of the 
filaments (controlled by the actin density param- 
eter Pact), filaments are only dilutely distributed 
and no network with intersections can consti- 
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(a) low pa 




cluster size 



Fig. 17. Cluster size distributions for different run- 
times. For 10000 time steps a stationary state is 
reached. 



(b) high Ud 



.• - •V'^ • 





(c) p = 0.008 




Fig. 16. Particle configurations (neighborhoods 
with \ — 2rp displayed) for system size L — 200 r^. 
Parameter examples where clustering is supressed. 
(a): Very low actin density pact — 0.04; no intersec- 
tions appear, thus no clusters can form, (b): High 
detachment rate ujd — 0.8: Particles are only rarely 
attached to filaments so that influence of the net- 
work is small and brownian dynamics prevail, (c): 
Low particles density; particles meet only rarely so 
that cluster cannot nucleate. 



tute. For small particle density events where par- 
ticles meet at intersections occur quite rarely so 
that cluster nuclei do not constitute. If the de- 
tachment rate approaches one or the attachment 
rate is low, certainly most particles are detached 
and perform a diffusive motion such that the in- 
fluence of the network vanishes. 



3.4.2 With branching 

In real actin networks, branching of filaments 
takes place quite frequently, resulting in a den- 
dritic network structure. If one is interested in 
the dynamics of vesicle transport on submem- 
branal actin networks, one has to consider this 
process as well. We checked CDs in a system with 
finite branching rate (here: branching probabil- 
ity uoi) = 0.0035) including the dependence of 
growth dynamics on the ARP2/3-density parp 
(cf. section [2^ . In this system, filament orien- 
tations are highly correlated. 

In fig. 21 we displayed cluster size distribu- 
tions for different particle densities p. 



Like in 

the more basic case of uncorrelated filaments, 
cj^ = 0, one observes an algebraic decay as well. 
This indicates that the scale free behaviour of 
the clusters is a robust feature of inhomogeneous 
transport networks and active particles exhibit- 
ing mutual steric interaction. 



3.4.3 No exclusion 

In order to distinguish the influence of exclusion 
and network structure, we switch off the exclu- 
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Fig. 18. Cluster size distributions in a diffusive 
system with a inhomogeneous active transport net- 
work. Default parameters are the ones given in ta- 
ble [4] while single parameters are varied: (a) system 
size L, (b) particle density p, (c) detachment rate 
cjcz, (d) actin density network density) pact- In 
a vast regime of parameter space the cluster distri- 
bution decays algebraically indicating the absence 
of a intrinsic size-scale. Only for small densities or 
extreme parameters clusters do not emerge. Default 
parameters: see table [4] 
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Fig. 19. Cluster size distributions in the inhomo- 
geneous network for different coarse graining scales. 
A is the radius of the neighborhood as defined in sec- 
tion |3]T] As well in the regime where large clusters 
emerge and in the non-clustering regime, the depen- 
dence on the coarse graining scale is weak indicating 
well separated clusters. 
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Fig. 20. Comparison of regular and inhomogeneous 
network displaying cluster size distributions for par- 
ticle density p = 0.04. While the CD decays expo- 
nentially in the regular network, its slope is alge- 
braic in the inhomogeneous one, demonstrating sig- 
nificant enhancement of clustering by the inhomoge- 
neous network structure. 
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Fig. 21. Cluster size distribution for branch rate 
(jjb > 0. Like in the system with non-correlated fil- 
aments, the CD exhibits algebraic decay for large 
clusters. 
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Fig. 22. Comparison of cluster size distributions 
with and without exclusion interaction. 



Fig. 23. Relative frequency distributions of particle 
densities in vicinity of intersection points. 



sion interaction. In figure 22 cluster size distri- 
butions of the inhomogeneous system with and 
without exclusion are compared for default pa- 
rameters. One observes that switching off exclu- 
sion suppresses formation of large clusters signif- 
icantly. 

In order to quantify the inhomogeneity of the 
particle distribution, we determined the average 
particle density in the proximity of intersection 
point^ Figure 23 displays the relative frequency 
of intersection points with local particles den- 
sity within density intervals of length 0.005. We 
see that while the density of particles is almost 
the same on all intersections in the regular net- 
works, there is a wide distribution of dilute as 
well as crowded intersection points in the inho- 
mogeneous network, marking a highly inhomo- 
geneous distribution of particles. 



3.5 Clustering mechanisms 

In this section we want to conjecture mechanisms 
involved in the formation of clusters in active 
transport systems. 

The most common clustering mechanism is 
probably nucleation of droplets in oversaturated 
van-der-Waals gases. If temperature or density 
is low enough, the attractive interactions of gas 
molecules favours aggregation of particles in clus- 
ters representing fluid phase droplets. A similiar 



^ Strictly speaking we virtually divided the system 
in boxes of edge length 3rp and determined its aver- 
age particle occupations considering only boxes that 
contain intersection points 



system was shown in section 3.2 to exhibit clus- 
ters on characteristic size-scales. 

Fig. [24] illustrates a mechanism that is si- 
miliar to clustering (jamming) at intersections 
in city street networks [17 . If two particles are 
on sites adjacent to an intersection site, they 
can block each other, forming an initial two- 
particle cluster. They can only leave the cluster 
by bypassing (which is inhibited for block rate 
6 = 1) or detachment followed by hopping on 
an adjacent non- filament site. The rate of these 
events corresponds to the flow out of the cluster 
J^^^. Particles attached to one of these filaments 
can run into this cluster from behind, contribut- 
ing to the flow into the cluster Jf^. If inflow is 
larger than outflow, a queue emerges similiar to 
high density regimes in TASEP-like driven lat- 
tice gases [9 . As long as inflow is larger than 
outflow, clusters grow, decreasing the number of 
free particles, thus Jf^ decreases. In contrast, if 
inflow is less than outflow, clusters shrink, re- 
leasing particles that contribute to an increas- 
ing free-particle density, which consequently in- 
creases Jf^. Finally, each cluster will reach a sta- 
tionary state with Jf^ = Jout^ providing a scale 
for the cluster sizes which is spatially invariant. 

In inhomogeneous networks intersection points 
also serve as nucleation seeds for clusters. How- 
ever, the density of intersection points may vary. 
If multiple filaments cross in a small area, clus- 
ters that extend over other intersection points 
constitute additional obstacles for particles mov- 
ing on other filaments. Hence clusters will ex- 
tend on these filaments (see fig. 25), where they 
might cover additional intersection points and 
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Fig. 24. Illustration of the mechanism that leads 
to the formation of clusters in a regular network. 
Particles blocking each other at intersections serve as 
a nuleation seed. Incoming particles are also blocked 
and increase the cluster size. 




Inflow: attachment+running into cluster 

Outflow: detachment, diffusing away 

Nucleation: 2 particles at crossing at the same time 



Fig. 25. Illustration of the mechanism that leads 
to the formation of clusters in a inhomogeneous net- 
work. If the density of intersections is high, existing 
clusters serve as additional obstacles for particles on 
other filaments, enhancing cluster formation. 



particles on many filaments can flow into that 
cluster. This way large clusters can emerge, while 
in regions with dilutely distributed intersection 
points only small clusters like in the regular net- 
work can constitute. 



4 Discussion 

We examined transport of hard core particles on 
regular and inhomogeneous networks embedded 
in a diffusive environment. The models consist 
of regions where particles perform diffusive mo- 
tion and of one dimensional stripes (filaments) 



to which particles can attach and perform di- 
rected motion. In most parameter regimes clus- 
ters are observed which are assumed to emerge 
due to blockage of particles at intersecting points 
of filaments. As a reference for clustering with- 
out transport networks we examined cluster fea- 
tures of a model with attractive interaction of 
particles. 

A detailed analysis of cluster size distribu- 
tions shows that the network structure as well 
as the particle-particle interaction change qual- 
itatively the form of the CD. On regular net- 
works one typically observes exponentially de- 
caying CDs at low densities, i.e. one can attribute 
a typical scale to the clusters. Algebraic CDs 
are only observed at large densities when single 
clusters merge to form large mesh-shaped cluster 
complexes. 

In contrast disordered networks, where fila- 
ment lengths and orientations are randomly dis- 
tributed, exhibit algebraic cluster size distribu- 
tions in a wide range of the parameter space, 
indicating that clusters on all size-scales exist. 
It is important to notice that the algebraic CDs 
are observed one order of magnitude below the 
respective densities in reference systems. In the 
disordered network, clusters are well separated 
exhibiting a compact structure. 

Switching off the hard core potential results 
in a strong suppression of large clusters, while 
the distribution of particles in the disordered 
network still remains inhomogeneous. Hence, in- 
homogeneity of the network structure as well as 
exclusion interaction strongly enhance the for- 
mation of macroscopic clusters. 

The analysis of cluster size distributions in 
the system of diffusive particles with attractive 
interaction shows clustering on a finite size-scale 
and a maximum in the cluster size distribution. 
While nucleation of clusters occurs on very small 
time-scales, the average size of clusters grows 
slowly. It is assumed that the stationary state is 
attained at very long times exhibiting only few 
or even a single large cluster including most par- 
ticles of the system. 

In summary, we have found that the micro- 
scopic particle dynamics as well as the network 
structure have significant influence on the quali- 
tative form of the cluster distribution. From our 
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point of view this observation is of great impor- 
tance for the analysis of biological systems since 
it is often not possible to identify the underly- 
ing microscopic mechanisms leading to exper- 
imentally observed aggregation. In these cases 
the analysis of the cluster distribution on larger 
scales may answer the question wether observed 
patterns are the result of active transport or of 
aggregation due to attractive interactions. 

We thank M. Schmitt for fruitful discussions and 
the German Science Foundation under grant number 
DFG GK 1276/1 for financial support. 
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